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Abstract 

We perform a detailed study of the phase diagram of the lattice Higgs 
SU(2) model with fixed Higgs field length. Consistently with previsions 
based on the Fradkin-Shenker theorem we find a first order transition line 
with an endpoint whose position we determined. The diagram also shows 
cross-over lines: the cross-over corresponding to the pure SU(2) bulk is 
also present at nonzero coupling with the Higgs field and merges with 
the one that continues the line of first order transition beyond the critical 
endpoint. At high temperature the first order line becomes a crossover, 
whose position moves by varying the temperature. 



1 Introduction 

In this work we present a detailed study of the phase diagram of the lattice Higgs 
SU(2) model with Higgs field in the fundamental representation and of fixed 
lengtrS The model in which Higgs length is allowed to change has received quite 
a lot of attention in the past for its possible phenomenological implications (see 
e -g- G±]) i so that its phase diagram is known with good precision. No systematic 
study exists instead for the case in which the Higgs length is frozen. However 
this last model is often used, mainly because of its computational simplicity in 
numerical simulations, as the prototype of a non-Abelian gauge theory coupled 
with matter in the fundamental representation. In particular some work has 
been done to study confinement in this model ([3] and [3]). In those works some 
properties of the phase diagram are usually taken for granted, like the existence 
of a line of first order transitions for ft > 2.3, but they have never really been 
tested in simulations. As a first step towards a complete understanding of this 
model, we thus started to systematically investigate its phase diagram, in order 
to obtain precise estimates of the location of its critical points. 
The action of lattice Higgs SU (2) model we adopt is 

S = P J2 I 1 - ^ReTriV(x) j - | Yl ^Wi^U^x + fi)</>(x + A)] (1) 



o v ' j 2 

where the first term is the standard Wilson action and the Higgs field <fi (which 
transforms in the fundamental representation) is rewritten as an SU (2) matrix 
(see e.g. [5], [S])- Since the action is linear in the fields at each point, standard 
hcatbath ( [7] , [5] ) and overrelaxation ( [5] ) algorithms can be used for the Monte 
Carlo update. If the Higgs length is allowed to change this is no longer true, since 



1 A first report of this study was presented at the LATTICE 2008 conference, £Q. 
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in this case a quartic term is needed, which has the form A{~Tr[<fy (x)cf)(x)] — l} 2 
and destroys linearity. 

In the limit /3 — > oo the theory reduces to an 0(4) non-linear sigma model, 
which is known to have a mean field phase transition for k w 0.6 (see e.g. |10|): 
in it was shown that this second order transition becomes a first order a la 
Coleman- Weinberg when the gauge field is introduced as a small perturbation (/3 
large). The authors of [TTJ, using the cluster expansion developed in [T2], were 
also able to prove the existence of a region of parameter space near the axis j3 — 
and for k — > oo where every local observable is analytic. This statement is often 
referred to as the Fradkin-Shenker (FS) theorem. It is important to notice that, 
in this context, an observable is defined as local when its support is contained 
in a compact set in the thermodynamic limit; observables not satisfying this 
requirement can have non-analytic behaviour (see e.g. [131 E]). This two results 
suggest a phase diagram like that shown in Fig. [TJ the analyticity region is 
indicated by AR and is limited by the dotted line, the thick line is the line 
of first order transitions and the two dots are its second order endpoints. As 
long as we consider the model at zero temperature, the e-expansion predicts the 
endpoints to be in the mean field universality class. 






Figure 1: Phase diagram of the Higgs SU{2) model as predicted in [TT] . 



A similar phase diagram was observed in the works [2J, [T3] for the lattice 
Higgs SU(2) model with fourth-order scalar coupling A < 0.5, while in the 
model considered here A — > oo. Because of the supposed triviality of the </> 4 
model in four dimensions, A is expected to be a marginally irrelevant parameter 
and therefore the phase diagram not to change qualitatively for A — > oo (see 
e.g. |16|); however it was observed since long time that the first order transition 
gets weaker as A is increased, so that the phase diagram of Fig. [JJhas not really 
been checked at large values of A. 

After the seminal work ref. [17, on very small lattices, in ref. [TB] the 
observation of a double peak structure was reported at j3 = 2.3 on a 16 4 lattice, 
which however was probably only a consequence of the poor statistics, since in a 
later study, ref. |19j . no double peak was found at (3 = 2.3. There it was stated 
that "the system exhibits a transient behaviour up to L — 24 along which the 
order of the transition cannot be discerned" . In this paper we produce the first 
clear evidence of the line of first order transition and we obtain an estimate of 
the endpoint position. 
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2 Numerical results 



The obvious observables to look at for this system are 

• the gauge-Higgs coupling, G = i(Tr[0 t (x)[/ AI (a; + (i)<f)(x + fl)]) 

• the plaquettc, P = |(TrP^> 

• the energy density E = 6/3 P + 4kG 

Besides these natural ones, we monitored also the following observables: 

• the Polyakov loop P L (x) = ±Tr [flto 1 U (t, S)] , P L = yiEsW)) 

• the Z2 monopoles density, M = 1 — -h-^2i C cr C i where c stand for the 
elementary cube and <r c = ] f signTrP M „ 

The Polyakov loop behaviour is used as an indicator of confinement. The study 
of the Z2 monopoles is motivated by the similarity of the first order transition 
with the bulk transition of the SU(2) pure gauge theory, which is driven by 
lattice artefacts such as the Z2 monopoles. Both these points will be discussed 
more accurately in the following. 

Data were analyzed by using the optimized histogram method ([20 ) and the 
statistical errors were estimated by using the moving block bootstrap method 
(see e.g. [21]). 

The presentation of the simulation results will be divided in several parts 

1. we will show that at /3 — 2.5 there is no signal of a phase transition and 
data are consistent with a smooth cross-over 

2. we will show that for j3 > 2.775 the scaling is consistent with a first order 
transition and we will present evidence of a double peak structure 

3. two independent estimates of the endpoint will be obtained 

4. we will give hints that the above transitions are not related to confinement 

5. we will investigate the relation between the line of first order transition, 
which becomes a smooth cross-over beyond the endpoint, and the pure 
SU(2) bulk transition 

6. finally we will present some exploratory results at T ^ 

All results in the first four parts have been obtained by fixing the value of j3 
and by looking for transitions in k. 



2.1 Cross-over region 

In Fig. [2] the maxima of the susceptibilities of G, P and M are plotted for 
various lattice sizes and (3 — 2.5. For lattices up to L ps 18 the maxima of 
susceptibilities are well described by a function of the form a + bL 4 , so that 
they seem to scale linearly with volume as expected for a first order transition. 
However on larger lattices all susceptibilities saturate and no singularity seems 
to develop in the thermodynamical limit. This means that the system has a 



3 




Figure 2: Top: G susceptibility peak heights at (3 — 2.5 (the continuous line is 
a fit to a + bx 4 ). Center: P susceptibility peak heights at /3 = 2.5. Bottom: M 
susceptibility peak heights at = 2.5. Clearly they all saturate at large L. 
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correlation length of order w 10 lattice spacing, so that the increase of the 
susceptibilities with volume, that in previous studies was interpreted as due to 
a first order transition, is just a signal that the lattices were too small. 

To look for transitions, besides susceptibilities we also check the behaviour 
of the energy Binder cumulant, which is defined as B = 1 — (E 4 ) / (3(E 2 ) 2 ). It 
can be shown (see e.g. [22 ) that near a transition B develops minima whose 
depth scales as 

B U = 3-3 (^W-'HO,^) (2) 

where E± = lim^o^fl), A = £J + — £L and e = 5(^4. + E-). In particular 
the thermodynamical limit of B\ m i n is less than 2/3 if and only if a latent heat 
is present. 

The scaling of the energy Binder cumulant minima at (3 — 2.5 is shown in 
Fig. [3J Also here two different behaviours are clearly visible: on small lattices 
there is scaling consistent with a first order transition. By increasing the volume 
B — > 2/3, indicating a smooth cross-over or a second order transition. 




v 1 

Figure 3: Binder cumulant minima of E at /3 = 2.5 (V = L 4 ). 

These results indicate that at ft = 2.5 there is no transition; this is in sharp 
contrast with all previous studies of this model, which concluded that for /3 > 2.3 
the system undergoes a first order transition. This wrong conclusion was based 
on the analysis of lattices of size up to 25 , which we have just shown to be too 
small. 

2.2 First order region 

At f3 = 2.775, f3 = 2.79 and /3 = 2.8 the scaling of the susceptibilities and 
of the energy Binder cumulant remains consistent with first order also for the 
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Figure 4: Scaling of the maxima of the susceptibilities and of the E Binder 
cumulant minima for f3 = 2.775 (V = L 4 and the larger V values correspond to 
lattice sizes L = 30, 35, 40, 45). 
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Figure 5: Histogram of the observable G for j3 = 2.8 on the two largest lattices. 



larger lattices, as shown for example in Fig. @]for j3 = 2.775. In this range of 
j3 values the transition gets stronger as j3 increases. However we know that the 
transition at j3 — > oo has to become of second order, so that going at (3 high 
enough the transition has to get weaker. We could not reach this regime in our 
simulations since when increasing the /3 value it is also necessary to use larger 
lattices. To have results free of spurious finite size transitions the lattice must 
be large enough for the corresponding pure SU(2) gauge theory to be confined 
and exponentially large lattices in /? are needed. 

Although all observables scale consistently with a first order transition, a 
clear signal of metastability was revealed only at j3 — 2.8, where the transition 
is stronger, and only in the two largest lattices, namely L = 45 and L = 50; 
the histograms for the observable G on these two lattices are shown in Fig. [5J 
where the formation of a double peak structure is visible. 

2.3 Endpoint 

Having three sets of data with first order scaling and knowing the universality 
class of the endpoint, we can estimate its position /3 C assuming that we are close 
enough to it. We can do this in two independent ways: 

1 . From the second form of equation [2] we know the dependence of B on the 
discontinuity A and for a mean-field transition we have A oc i 1 / 2 , where t 
is the reduced temperature, t = (T — T c )/T c . In our case, near the critical 



7 




8 



point, we can use t oc (/3 — /3 C ) and therefore equation [2] can be written a^| 

lim B L \ mm = \~ 7l (/3 - p e ) + 0(/3 - /3 C ) 2 (3) 

L— too o 

2. For a first order transition the susceptibilities scale for large volume as 

X = V((0 2 ) - (O) 2 ) « const + VA 2 (4) 

so that fitting the x maxima to the relation \ — const + VA 2 and using 
again A oc t 1 / 2 , we obtain 

A 2 = 72 (,3 - Pc) + 0(f3 - /3 C ) 2 (5) 

The fit in the planes (j3,B) and (/?, A 2 ) are shown in Fig. [5] (we use the 
G susceptibility); in both cases the fit is very good and the estimates for the 
critical point position are 

= 2.7266(16) f3 {2) = 2.7299(36) (6) 

which agree with each other within errors. 

2.4 Relation to confinement 

The next question is to understand if the two "phases" found at large /3 have 
different confinement properties. We recall that they cannot be considered as 
different thermodynamical phases as a consequence of FS theorem. 

In this model Wilson loops never obey the area law because of the presence 
of the Higgs field, which destroys the center symmetry of the pure gauge theory; 
as a consequence the Polyakov loop is always non-zero (a direct check of this 
is shown in Fig|5]) and cannot be an order parameter. Nevertheless Polyakov 
loop is commonly used as a confinement tracker also in theories where it is not 
an order parameter, e.g. in QCD, since it abruptly jumps at the deconfinement 
transition. For the lattice SU(2) Higgs model the Polyakov loop doesn't seem to 
be influenced in any way by the transition: for small lattices it slightly increases, 
however for larger ones, the transition gets stronger but Pl gets flatter in the 
transition region, as can be seen in Fig. [5](the transition is at k = 0.704675(30)). 
Also Polyakov loop correlators, measured by using the multilevel algorithm of 
ref. |24) . do not show any significant change across the transition, as shown in 
the bottom of Fig. [8] 

An alternative possibility is to use the order parameter introduced in [25], 
which we will denote by Ofm- This is defined by the limit for R — » oo of 
the quantity Ofm(R), which is constructed as the square of the mean value 
of a staple-shaped parallel transport connecting two Higgs fields divided by a 
Wilson loop; the size of the staple is related to that of the Wilson loop as shown 
in Fig. where Higgs fields are represented by dots. To our knowledge Ofm 

2 A subtlety has to be considered here: since in this model the two relevant operators at 
the mean field critical point are not related to any symmetry of the system, there is in general 
mixing between the magnetic and thermal relevant operators (see e.g. |23l ). However, since 
the A's are measured along the coexistence curve, near the critical endpoint the mixing is 
negligible and the relation between A and t is as usual A oc t@ (here "j3" is the ft critical 
exponent). 
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Figure 7: Ofm(R) definition. 



has never been measured in a Monte Carlo simulation, so we have no previous 
results to compare with; its strong coupling limit was computed in [26] . Ofm (R) 
measures the overlap between a Higgs-Higgs dipole of linear dimension R and the 
vacuum; if asymptotic colored states exist they are orthogonal to the vacuum, 
so Ofm{R) —> for R — >• oo, otherwise we must have Ofm(R) — > ce > 0. 
The results obtained by measuring Ofm{R) on a 40 4 lattice at ft = 2.775 for k 
slightly above and slightly below the transition points are shown in Fig. |H1 where 
the two lines are fits to the functiorjf] f(x) = a + b-^ (symmetrized because of 
periodic boundary condition). 

We have for the parameter a, c the estimates: 

a(o.7046) = 4.9(1) • 1(T 5 c (0 . 70 46) = 2.98(5) (7) 
a ( o.7048) = 10.5(5) • 1(T 5 c (0 . 7048) = 3.30(23) (8) 

so on both sides of the transition charge is screened and according to the inter- 
pretation of [25] there is color confinement. 

2.5 The pure gauge SU(2) bulk 

Another interesting point to study is what relation (if any) the first order has 
with the bulk transition of the SU(2) pure gauge (we call this "transition" for 
the sake of simplicity, although it is only a rapid but smooth cross-over). An 
hint on the existence of such a relation can be obtained by noting that the value 
ft ps 2.3 previously thought to be the first order line endpoint position is the 
value of the bulk SU (2) transition. 

Since some test simulations indicated that the position of the SU(2) bulk 
is very stable for small values of k (so that in the plane {ft, k) this crossover 
follows a line perpendicular to the ft axis), it seems more convenient to hold k 
fixed and vary the ft value; the results obtained on a 30 4 lattice for the plaquette 
susceptibilities for several k values are shown in Fig. 1101 

From these data a structure of cross-over lines emerges as shown in Fig. 1111 
For k < 0.6 the bump in ft of the plaquette susceptibility which corresponds to 
the SU(2) bulk transition is independent on n both in position and in shape 
(vertical dotted line in Fig. [TT|) . For n larger than 0.6 two peaks in ft appear, the 
bulk one and the first order transition studied above. The latter peak persists 
also when k is greater than the critical value corresponding to the endpoint 

3 This ansatz is motivated by the observation that for R big enough such that the Wilson 
loop scales with perimeter, the exponentials in numerator and denominator are the same, 
leaving a power law. 
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Figure 8: Top: Polyakov loop for /? = 2.775, «; = 0.70464 and various L, 
the line is a fit to a + fecxp(— cL) and the result for the asymptotic value is 
a = 1.7(1) • 10~ 4 . Center): Polyakov loop for (3 = 2.775 and various L. Bottom: 
Polyakov loop correlators measured on a 40 4 lattice for (3 = 2.775 and k slightly 
below (0.7046) and above (0.7048) the critical value. 
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Figure 9: Results for Ofm(R) at f3 = 2.775 on a 40 4 lattices slightly above and 
below the transition point. 



andthcre is no real transition, although it is smaller that the SU (2) bulk signal 
(see Fig. [TU1 k — 0.75). By increasing the k value the cross-over remnant of 
the first order transition moves towards smaller j3 values, until it intersects the 
SU(2) bulk in the point indicated by A in Fig. [TTJ In the neighborhood of 
the point A the "first order continuation" peak gets stronger and the SU(2) 
bulk disappears. This increase could have been misinterpreted as the first order 
transition line endpoint in previous works. A direct check to ensure that in this 
region there is no transition is shown in Fig. 1121 For still larger k values only 
one maximum is present in susceptibilities, which gets weaker as n — » oo. 

This interplay between the first order transition line and the bulk SU (2) 
suggests that the first order transition could be in some way related to the same 
lattice artifacts that drive the bulk SU (2) transition, namely the monopoles. 
The density of Z2 monopoles seems indeed to have a jump across the transition 
(see Fig. gj|. 

2.6 Finite temperature 

Motivated by the last observation we try to investigate if the first order transi- 
tion line can itself be thought as a lattice artefact, like the pure gauge 5*0(3) 
first order bulk transition. To answer this question the simplest method is to 
consider the system at finite temperature: bulk transitions are insensitive to the 
temporal extension L t of the lattice, while physical transitions scale by varying 



For each L t value we simulate the system with several spatial extents L s , in 
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Figure 11: Phase diagram of the Higgs SU(2) with rapid crossovers indicated 
by dotted lines. 

order to perform a finite size analysis. Some of the results obtained are shown in 
Fig. [T3J It is clear that the position of the maxima of the G susceptibility moves 
by varying L t , so that the line of first order transition seen at zero temperature 
cannot be a lattice artifact. Moreover for all the values of L t considered the 
maxima saturate by increasing L s , indicating the absence of a phase transition. 
We expect that for sufficiently small temperature (i.e. large enough L t ) the first 
order transition is restored, however to observe this behaviour numerically we 
should use lattices with L t > 20 and L s 3> L t , typically L s > 3L t . Instead we 
have to restrict ourselves to L\L t <4x 10 6 because of computer capability. In 
any case we see the cross-over gets stronger by increasing the lattice temporal 
extension. 

3 Conclusions 

This paper is a study of the phase diagram of the SU(2) Higgs gauge theory with 
the Higgs in the fundamental representation and fixed length. We find that to 
investigate the phase structure of this system it is necessary to use much larger 
lattices than the ones adopted in the past due to the large correlation length. 

We present the first clear evidence of the presence of a line of first order 
transition and we estimate its endpoint position to be at p c — 2.7266(16), thus 
confirming an expectation based on the Fradkin-Shenker theorem. 

We give indications that this transition is not a deconfincment transition. 

We attempt an exploration of the system at finite temperature; the first 
order transition becames then a cross-over whose position moves by varying the 
temperature. 
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Figure 12: Plaquette susceptibility for k = 0.85. 



A Data 

In this appendix we give some details of our numerical results. The notation is 
as follows 

• K pc is the location in k of the G susceptibility maximum 

• xo is the peak value of the susceptibility of the O observable 

• B is the value at minimum of the energy Binder cumulant 



= 2.725 


L 


ftpc 


Xg 


Xp 


XM 


XE 


2/3-B 


20 


0.7090(2) 


0.615(18) 


0.005575(88) 


0.1147(19) 


10.85(30) 


6.36(18)e-07 


25 


0.7089(2) 


0.692(16) 


0.005831(71) 


0.1194(15) 


11.89(26) 


2.839(62)e-07 


30 


0.7088(1) 


0.858(20) 


0.006485(78) 


0.1319(15) 


14.51(32) 


1.679(38)e-07 


35 


0.70873(5) 


1.018(24) 


0.007082(89) 


0.1431(17) 


16.64(37) 


1.039(23)e-07 


40 


0.70866(5) 


1.022(24) 


0.00704(12) 


0.1426(18) 


16.85(37) 


6.17(13)e-08 


45 


0.70886(5) 


0.994(31) 


0.00725(23) 


0.1469(25) 


16.25(47) 


3.71(10)6-08 
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/3 = 2.775 


L 




XG 


Xp 


Xm 


XE 


2/3-5 


25 


0.7049(1) 


0.5740(65) 


0.005459(31) 


0.10807(62) 


9.460(98) 


2.156(22)e-07 


27 


0.70480(8) 


0.6211(82) 


0.005528(40) 


0.10884(74) 


10.16(12) 


1.702(21)6-07 


30 


0.70476(5) 


0.734(13) 


0.005769(52) 


0.1135(10) 


11.82(19) 


1.299(21)6-07 


35 


0.70466(5) 


0.812(16) 


0.006044(64) 


0.1182(12) 


12.88(24) 


7.64(14)e-08 


40 


0.70469(3) 


0.956(17) 


0.006573(77) 


0.1285(12) 


14.66(25) 


5.099(88)e-08 


45 


0.704675(30) 


1.114(24) 


0.00698(10) 


0.1359(16) 


16.72(34) 


3.631(75)c-08 



P = 2.79 


L 


ftpc 


XG 


Xp 


XM 


XE 


2/3-5 


25 
30 
40 
45 


0.70380(10) 
0.70360(5) 
0.70356(3) 
0.703525(30) 


0.5512(97) 
0.646(12) 
0.912(19) 
1.140(23) 


0.005301(48) 
0.005481(49) 
0.006187(78) 
0.006901(97) 


0.10362(93) 
0.10632(95) 
0.1200(13) 
0.1318(15) 


9.05(14) 
10.39(17) 
14.01(27) 
17.01(33) 


2.030(32)c-07 
1.125(19)e-07 
4.799(93)e-08 
3.639(70)c-08 



/3 = 2.8 


L 




XG 


Xp 


XM 


XE 


2/3-B 


25 


0.7032(2) 


0.683(22) 


0.005080(80) 


0.0978(15) 


11.25(34) 


2.500(75)e-07 


30 


0.70298(8) 


0.7212(85) 


0.005279(38) 


0.10181(60) 


11.63(13) 


1.246(13)e-07 


35 


0.70286(5) 


0.760(16) 


0.005685(62) 


0.1093(11) 


11.64(22) 


6.74(13)e-08 


40 


0.70282(4) 


0.946(22) 


0.006239(75) 


0.1192(13) 


14.36(31) 


4.87(ll)e-08 


45 


0.70279(3) 


1.219(27) 


0.007033(89) 


0.1339(16) 


17.75(37) 


3.759(78)e-08 
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